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Abstract 

The classical transverse field Ising spin- glass model with short-range inter- 
actions is investigated beyond the mean- field approximation for a real d- 
dimensional lattice. We use an appropriate nontrivial modification of the 
Bethe- Peierls method recently formulated for the Ising spin- glass. The zero- 
temperature critical value of the transverse field and the linear susceptibility 
in the paramagnetic phase are obtained analytically as functions of dimen- 
sionality d. The phase diagram is also calculated numerically for different 
values of d. In the limit d — > oo, known mean- field results are consistently 

reproduced. 
75.10.Nr, 75.50.Lk 
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The study of glasses is today one of the most relevant and actual problem in condensed 
matter physics. Originally, the basic idea was to start from spin glass (SG) models and 
to extract as much as it was possible at a mean- field- approximation (MFA) level |3]-fIJ. 
However, there are recent studies |5|-|8| which indicate difficulties to extend the MFA scenario 
to realistic spin glasses with short- range interaction and decide " a priori" which properties 
survive and which must be appropriately modified. Renormalization group treatments [PHTT 



for classical and quantum spin glasses and phenomenological studies |L2[ do not seem to 
suggest a clear picture. 



Quite recently, in an interesting paper ||13|| , an approach beyond the MFA has been 



achieved for an d- dimensional Ising SG model with short- range interactions on a real 
lattice using an extension of the Bethe- Peierls approximation (BPA) |14| to the spin glass 
problem via the replica trick. This approach seems to be very promising to estabilish a direct 
contact with the results obtained by different authors for the infinite- ranged version and to 
controll possible deviations for short- ranged glasses from the well acquired MFA scenario. 
Of course, additional applications to more complex glassy systems and improvements are 
necessary for understanding something more about the role played by the glassy fluctuations 
around the MFA solution for finite dimensionalities. 

In this paper we explore the glassy properties of the d- dimensional classical transverse 
field Ising SG model [|H)Hl7|1 with short-range interactions using an appropriate nontrivial 
modifications of the BPA formulated originally in Ref. [|13[ for the Ising SG. The model here 



considered has received recent attention because it is a relatively simple SG model which 
reflects some properties of the quantum counterpart [|] and it is a specific example of classical 
two- vector anisotropic SG. So, the classical limit of the usual more complex quantum Ising 
SG model in a transverse field with realistic exchange interactions may be useful for taking 
contact Hn||r?| with low- temperature properties of the so- called "proton glasses" |15[, as 



the compounds Rbi.^ (NhLi)^ H2PO4, and with the most recent experimental magnetic data 



for the dipolar glass Li H00.167Y0.822F4 [P~8||T9 



As concerning the quantum realistic SG model, a lot of results have been obtained only 



for d = 1 pO| , |2lH and for infinite- ranged interactions (d = oo) For the classical coun- 

terpart, only results with infinite- ranged interaction have been derived |IB|JI7|| , in particular 
at T = 0, and even in the "simple" MFA limit the full phase diagram has not yet been cal- 
culated. In any case, there are rather little studies about short- ranged glassy models. So, it 
appears quite relevant that the BPA allows us to describe some nontrivial glassy properties 
of the model (in the paramagnetic phase) for d > 1. In particular, explicit analytical results 
are obtained at T = and the phase diagram in the (temperature, transverse field)- plane 
is derived numerically for an arbitrary dimensionality. 

The classical transverse field SG model here considered is described by the Hamiltonian 



IM3: 
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i=l 



with (Sf ) 2 + (S^) 2 = l.Here T and h are transverse and longitudinal fields, respectively, 
the couplings Ji/s are independent random variables assuming values ±J with the equal 
probability. In ([[]) Y,uj) ' ' ' denotes a sum over nearest- neighbours pairs of N sites on a 
hypercubic d- dimensional lattice. Using the replica trick, all is reduced to determine the 
"quenched average": 



Trexp [-(3^2 H a 



a=l 



(2) 



where H a is the a-th replica of the Hamiltonian ([I]) and (3 = 1/T with the Boltzmann 
constant k-Q = 1. Working directly on the real lattice, the basic idea of the BPA for spin 
glasses |TB[ is to take into account the correct interactions inside replicated clusters (cl), 
constitued by a central spin So and its 2d nearest- neighbours {Si; i = 1, • • • ,2c?}, and to 
describe the interactions of the cluster borders with the remnant (rm) of the system by 
means of effective couplings among replicas to be determined self - consistently. With this 
in mind, Eq. (^j) for the Bethe- Peierls ansatz can be formally revritten as [13[ 

Z n = Tr {Scl} exp (-0 f) H<A Tr {Srm} exp (-/3 £ H t 
= K{T,Y,h) Tr {Sci} [ex P (-/3W n )] av 
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where 

2d 2d 

Ha — _ r ^ Sl a - ^2 Joi^Oa^ia (4) 
k=0 i=l 

and H^ m ^ denotes replicated Hamiltonians of the cluster and remnant of the system inter- 
acting with cluster borders, respectively, and K (T, T, h) is a multiplicative constant inde- 
pendent on lateral spins, 

1 A 2d n jqz 

and 

n 2d BJ^ n ^ 

7~Ln = - 22z2 JoiSoaSia 7T E E ^aa' ^fa^ta' (6) 
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with A aQ / = /i aa / for a ^ a' and Aq, q = /i which are parameters to be determined via 



appropriate self- consistent equations. Here we have used the relation S% a = ±y 1 — (S^ a ) 2 
(k = 0,1,- ■■.2d). Of course, if a transition from a paramagnetic phase to a SG one is 
assumed to exist, one expects n aa > = in the paramagnetic phase. 

At this stage, the self- consistent equations which determine the effective couplings fi a>a i 
and /i as n — > 0, are 

(S* a S? al ) = (S* 0a S* al ) with i = 1, • • ■ , 2d , (7) 

where 

= Tr [exp (-g^Q ■ ■ -] av 
^ ' Tr [exp (-3H n )l v ' 1 J 

It is easy to check that, for h = 0, due to the inversion symmetry Sf a — > — S^a and symmetry 
of the probability distribution for J^- Eq. (|7|) with a = a' can be reduced to the following 
one: 

Xi = Xo for % = 1, • • - ,2d , (9) 
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where 



X* 



d(S ka ) 
dh 



=0 with k = 0, 1, • • • , 2d 



(10) 



denotes the local susceptibility. 

We are now in the position to obtain the explicit equations for fi which will be used 
for obtaining also the phase diagram of the model. Since it is expected that \i a ^ a i — > 
approaching the spin- glass transition from below, for n — > 0, one obtains at h = 

2d 



(Shasta') — P^Vaa' ^ I ( $i ) o] av + O 



i=l 



(k = 0, 2d) 



(11) 



where 



with 
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fc=0 



(12) 
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In Eq. fljjD denotes the normalization factor. The term —h Y^Lq SI must be added 



to the right hand side of (0) when it is necessary. So, at h — due to the translational 
symmetry for the sample averaged system and assuming \i aQ i = at and above the glassy 
transition line (to be determined), the self- consistent equation (0) for a^a' and a = a' 
reduces Respectively, to: 



((St) 2 )' + (2d-l) (S*S*)l 



2d 



(SoSj) 



(14) 



2d) 



and 



((S7) 2 )c 



((S*) 2 ) (i = l,...,2d), 



(15) 



where z ^ j denote arbitrary lateral sites of the cluster with the central spin So- 



By solving Eqs (0)- (|15D , it is possible to obtain the phase diagram of our model in 
the (T, T) plane. Explicit results can be derived analytically only at T = 0. As T — ► 0, 
introducing /Z = /3fi which is finite, with the help of Eq. ([12]) choosing in Eq. i = 1 one 
obtains: 



1 r- Jr 2 -4(2c/- 1) J 2 



(16) 



J 2 2 

and hence taking into account Eq. @ one finds, in the paramagnetic phase at T = 0, linear 
susceptibility: 



X = 2 



r + ,/r 2 - 4 (2d - 1) j 2 



(17) 



T + v /r 2 -4(2rf-l) J 2 j -4J 2 
Now, we calculate the critical value r c at T = of the transverse field using Eqs (14) and 



(|16|). For i = 1 and j ^ 1 , with some algebra we rewrite Eq. (jT4|) , as T — > 0, in the following 
form: 



r c - /i c J 2 ) 4 - 2c/ j 2 (r c - iz c j 2 ) 2 + (2d - 1) j 4 = o 



(18) 



with /x c J 2 = | T c — wT 2 — 4 (2d — 1) J 2 . From this equation one easily obtains: 



T c = 2 (2d- 1) 1/2 J . 



(19) 



As we see, % is positive and has physical meaning only for T > T c . This suggests that the 
expression (|TTj) for x is related only to the paramagnetic phase. 
For (r - r c ) /r c < 1 Eq.(0) yields: 



X « < 



(2d-l) 1 / 2 




2J(d-l) 




1 ( 









l + O (I* 



1/2 



for d ^ 1 
for fi = 1 



(20) 



. As an "a posteriori" justification of the correctness of the glassy BPA (|3]), it is easy to 
check analytically that, using the rescaling J — » J j \[2d one finds /Z = 2cZx and we get at 
T = for d ^ oo; x = 2J^ [r - (r 2 - T 2 ) 1/2 ] with T c = 2J + O {d' 1 ' 2 ). These results 
reproduce exactly those obtained at T = for the same SG model but with infinite- ranged 
interactions |l(J. This partial result supports the validity of the BPA for SG's. 
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The situation for d = 1 with a divergence of the linear susceptibility at T = as 
r — > can be simply explained. With the dichotomic probability distribution of one - 
dimensional nearest- neighbours couplings J (i, i + 1) = Jj, ( Jj = ± J with J > 0) after the 
gauge transformation of spin variables 



the sytem can be reduced to the uniform ferro magnet in an external transverse field T. 
Therefore it is naturally to expect that at T = T c the ferromagnetic phase transition with a 
divergent linear susceptibility occurs. Indeed a more detailed analysis of the one dimensional 
case shows that at T = the linear susceptibility x can be calculated exactly for the 
paramagnetic phase. The divergence of x is the same as that obtained within the BPA for 
the one- dimensional system. 

For arbitrary T and T, from Eqs (|14]) and (|15|) one can calculate numerically equilibrium 
properties ofour model in the paramagnetic phase. In particular in Fig.|l|, the phase diagram 
in the (T, T)- plane for different d is shown. We have conveniently scaled variables T and T 
for reproducing results at very high dimensionality. In Fig. ^| a variation of rescaled critical 
temperature with a dimension at T = is plotted. 

In conclusion, we have studied some relevant aspects of the classical transverse field short- 
ranged Ising SG in the paramagnetic phase for arbitrary dimensionality d. We expect that 
our results may be also useful for explaining some properties of the quantum counterpart of 
the model here considered. However, some questions remain to be explained. For example, 
on the basis of the general self- consistent equation (0) it is interesting to find solutions with 
fi aa i in order to see if the BPA is able to describe correctly our model in the SG phase 
at arbitrary dimensionalities. Within present calculations working for paramagnetic phase 
this is practically impossible, since the complicated integral ([12]) has been reduced at T = 
to the asymptotic form being Gaussian like one. Such an asymptoptic form is insufficient 
when parameters included even in the replica symmetric form. Therefore further 

works will be necessary to elucidate these problems. 




(21) 
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FIGURES 

FIG. 1. Phase diagram of the classical transverse Ising spin- glass with short- range interaction 
within the Bethe- Peierls approximation for spatial dimensions d = 2, 3, 4, 5 and 6. The temperature 
T and transverse field T are rescaled by the factor (2d) 1 . The larger the dimension, the higher 
the corresponding line. Here J = 1. 

FIG. 2. Rescaling critical temperature T c (2d)~ 1 ^ 2 for T = versus the dimension d. Here J = 1. 
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